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Using exact computations we study the classical hard-core monomer-dimer models on m x n plane 
lattice strips with free boundaries. For an arbitrary number v of monomers (or vacancies), we found 
a logarithmic correction term in the finite-size correction of the free energy. The coefficient of the 
logarithmic correction term depends on the number of monomers present (v) and the parity of the 
width n of the lattice strip: the coefficient equals to v when n is odd, and v/2 when n is even. The 
results are generalizations of the previous results for a single monomer in an otherwise fully packed 
lattice of dimers. 
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The monomer-dimer problem is a classical model in 
statistical physics. In the model the diatomic molecules 
are modeled as rigid dimers which occupies two adja- 
cent sites in a regular lattice and no lattice site is cov- 
ered by more than one dimer. The lattice sites that are 
not covered by the dimers are regarded as occupied by 
monomers. A central problem of the model is to enu- 
merate the dimer configurations on the lattice. A break- 
through came in 1961 when an elegant exact closed-form 
solution was found for a special case of the model in plane 
lattices, namely when the lattice is completely covered 
by dimers (the close-packed dimer problem, or dimer- 
covering problem) The method used to achieve this 
solution, however, cannot be extended to tackle the more 
general monomer-dimer problem where there are vacan- 
cies in the lattice, and the problem remains unsolved 
notwithstanding years of efforts. Rigorous results ex- 
ists only for series expansion in low dimer density 0, 
lower bounds on free energy |3j , monomer-monomer cor- 
relation function of two monomers in a lattice otherwise 
packed with dimers Q, and locations of zeros of parti- 
tion functions [j| . Some approximate methods have been 
proposed One recent advance is an analytic solution 
to the special case of the problem where there is a single 
vacancy at certain specific sites on the boundary of the 
lattice pj . 

There is renewed interest in monomer-dimer problem 
recently. The interest comes from different directions. 
Besides the intrinsic interest of the problem itself and 
its close relation to the well-studied Ising model [g, the 
model also acts as the classical limit of the recently intro- 
duced quantum dimer model [jj, which has been investi- 
gated intensively as the central model in modern theories 
of strongly correlated quantum matter. 

The problem also attracts attention in the field of 
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computational complexity. It has been shown that two- 
dimensional monomer-dimer problem is computationally 
intractable and belongs to the "#P-complete" class |l(l . 
#P-complete class plays the same role for counting prob- 
lems (such as counting dimer configurations discussed 
here) as the more familiar NP- complete class for the de- 
cision problems (such as the well-known traveling sales- 
man problem). The #P-complete problems are at least 
as hard as the NP-complete problems If any prob- 
lem in the #P-complete class is found to be solvable, 
every problem in #P class is solvable. Currently it is not 
clear whether there exists any such solution to the #P- 
complete or NP-complete class problems, and "P verse 
NP" problem is the perhaps the major outstanding prob- 
lem in theoretical computer science. 

In this Letter we report the logarithmic correction in 
the free energy of the monomer-dimer model with arbi- 
trary number of monomers on plane lattice strips of size 
ttl x 7i, where the width n is fixed, and address the finding 
in the context of universality, scaling, and exact finite- 
size corrections in critical systems. These topics have 
been studied intensively in recent decades. As one of the 
few non-trivial exactly solved statistical lattice models, 
the dimer model attracts much attention recently to test 
the predictions of conformal field theory [12( and finite- 
size scaling ^jjj]. Unlike other models, the (close-packed) 
dimer model on plane lattice shows certain peculiarity 
with respect to the predictions of the theories, both for 
the central charge |14L fl5| and for the logarithmic correc- 
tions |l4j . The free energy of a two-dimensional system 
at criticality is assumed to have the form of [lj, 

F m ,n = mnf b + (m + n)f s +Chxm + D + o(l). (1) 

The first two terms, ft, and f s , are the "bulk" and "sur- 
face" terms, respectively. These two terms are not uni- 
versal and depend on the details of the model. The other 
two coefficients, C and D, are supposed to be univer- 
sal, depending only on the shape of the system and the 
boundary conditions. The dimensionless coefficient D is 
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related to the central charge of the system and the bound- 
ary conditions in the transversal direction |12| . The co- 
efficient of the logarithmic term C depends on the geom- 
etry of the system. For systems with corners, as for the 
lattices with free boundaries discussed here, C = Y"! - m , 
with some universal contribution Ui from each corner . 

From the exact solution of the close-packed dimer 
problem, the asymptotic expression of the free energy 
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a mn /2 was found to be 



F m , n = mnf b +(m+n)f s +D n [(m+l)/(n+l)]+o{n 2+s ) 

(2) 

where 5 > and D n depends on the geometry of the 
lattice and the parity of n. Two issues arise from this 
asymptotic expression. The first issue is related to the 
central charge, because the expansion depends on the 
parity of n: 
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when n is odd. Here G is the Catalan's constant. At the 
first look, the coefficients of n~ 2 , which are related to the 
central charge of the conformal field theory, are different 
for even and odd n, leading to c = 1 for the even n and 
c = — 2 for odd n. Recently this issue has been discussed 
using logarithmic conformal field theory 0, ^J. It is 
pointed out that the coefficients of n~ 2 are actually re- 
lated to the effective central charge c c ff = c — 24/i m ; n , 
where h m i n is the conformal weight of the ground state 
and is a boundary dependent quantity. By using a bi- 
jection of the dimer coverings with a spanning tree and 
Abelian sandpile model, it is found that changing the par- 
ity of n results in the change of the boundary conditions 
of the mapped spanning trees, leading to different c e s- 
The central charge itself, however, remains unchanged 
c=-2 0. 

The second issue, which was pointed out previously 
|14| , is the absence of corner contribution term predicted 
by the theory (Eq. in the expansion of Eq. (J2J 
of the close-packed dimer model. This issue, however, 
changes its nature when we look at the general monomer- 
dimer model. Recently we found that when there is one 
monomer in the lattices, there is indeed a logarithmic 
term in the finite-size correction of the free energy |18| . 
In the following we generalize the results to plane lattices 
with arbitrary number of monomers, and found that the 
coefficients of the logarithmic correction term depend not 
only on the number of monomers, but also on the parity 
of the width of the lattice strip n. 

The configurational grand canonical partition function 
of the system is 



where is the number of ways to arrange k dimers on 
the mxn plane lattice with free boundaries. We are inter- 
ested in the coefficients of fixed number v of monomers, 
that is, a,(mn~v) /2i when the width of the lattice strip n 
is fixed and the length of the strip m changes. When the 
width of the lattice strip n is even, v can only take even 
numbers; when n is odd, v can take either even or odd 
values based on the parity of m. 

The full partition function have been calculated ex- 
actly for n from 1 to 16, using the method discussed in 
Refs. 0,0. The coefficients a( m „_ 1 ,)/2 of the partition 
functions are extracted to fit the following function: 
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where b = 1 when n is odd, and 6=1/2 when n is 
even. The reason to choose ln(m + 1) instead of ln(m) is 
discussed in Ref. [l^, an d the choice will not affect the 
results discussed below. 

For the close-packed dimer problem (where v = and 
mn is even), the exact solution gives Co in Eq. as []J 
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For a given n, this exact expression of cq actually holds 
for any finite number v of monomers. This was confirmed 
by fitting the data with cq as a free parameters. In the 
following we use cg(n) given in Eq. Q for c$. We fit the 
data extracted from the full partition functions to Eq. 
© for v = 0, . . . , 12 for odd n, and for v = 0, 2, 4, ... , 24 
for even n. The fittin g re sults for I are shown in Tables 
U and |nj respectively |20j . As we did previously 0] , in 
these fittings, only m > mo = 100 are used. The curves 
of fitting for two values of n, n = 11 and n — 12, are 
shown in Figure ^ and Figure |3 respectively. The fitting 
results for all the parameters when n = 11 and n = 12 are 
shown in Tables II I II and IIVI respectively. These results 
lead clearly to the conclusion that for a fixed number of 
monomers v, there is a logarithmic correction term in 
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of this term not only depends on v, but also on the parity 
of n, the width of the lattice strip: when n is odd, the 
coefficient equals to v; when n is even, the coefficient is 
v/2. The lack of the logarithmic correction term found 
in the close-packed dimer problem (v = 0) Q, Q 
and the logarithmic correction term found in the odd- 
by-odd lattice with one single vacancy (v = 1) ^| are 
just special cases of this general result. 

The close-packed dimer model, which is a special case 
of the general monomer-dimer model discussed here, al- 
ready demonstrates intriguing connections with the the- 
ory of university and finite scaling 0, ^| . It would be 
interesting to see how the results reported here incorpo- 
rate into the bigger picture. 
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TABLE I: Coefficient of the logarithmic term I from fitting (mn)" 1 ma( mn _„)/2 to Eq. [3]for odd values of n, with 6=1. Only 
data with m > mo = 100 are used in the fitting. For n = 1, m < 2000; for n = 5, . . . 13, m < 500; for n = 15, m < 200. 
Numbers in square brackets denote powers of 10. 
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6.98911-10] 


7.10235[-10] 
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1.56222[-9] 


3.76785[-9] 


3.42299[-9] 
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1 
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0.999999 


0.999998 


0.999997 


0.999993 


0.999968 


3 
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0.999998 


0.999991 


0.999975 


0.999945 


0.999897 


0.999781 


0.999088 
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0.999988 


0.999945 


0.999851 


0.99969 


0.999449 


0.998907 


0.996177 
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0.999962 


0.999835 


0.999576 


0.999163 


0.998585 


0.99736 


0.991937 
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0.999912 


0.999632 


0.999108 


0.998322 


0.997282 


0.995194 


0.986896 
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0.999836 


0.999351 


0.998494 


0.997273 


0.995726 


0.99273 


0.981618 
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0.999721 


0.998951 


0.997671 


0.995927 
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0.98982 
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0.999581 


0.998493 


0.996767 


0.9945 


0.991812 


0.986905 


0.970599 


10 
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0.99939 


0.997902 


0.995649 


0.99279 


0.989492 


0.983621 


0.965005 


11 
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0.999181 


0.997289 


0.994527 


0.991117 


0.987267 


0.980527 


0.959804 


12 


0.999999 


0.998909 


0.996529 


0.99318 


0.989159 


0.984713 


0.977081 


0.954453 



TABLE II: Coefficient of the logarithmic term / from fitting (mn)~ lna( mn _„)/2 to Eq. EJfor even values of n, with b — 1/2. 
Only data with m > mo = 100 are used in the fitting. For n — 2, m < 2000; for n = 4, m < 1000, for n = 6, . . . 14, m < 500; 
for n = 16, m < 130. Numbers in square brackets denote powers of 10. 
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1.216431-11] 


1.3329[-10] 


5.89405[-10] 
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0.999994 
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0.999739 
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0.993311 


0.98074 


10 
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0.999997 


0.999907 


0.999639 


0.999043 


0.998001 


0.987848 


0.977344 


12 


1 


0.999982 


0.999718 


0.999103 


0.997943 


0.996201 


0.983945 


0.981147 


14 


1 


0.999951 


0.9994 


0.998314 


0.996526 


0.994189 


0.982769 


0.991167 


16 
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0.999897 


0.998938 


0.997305 


0.994942 


0.99226 


0.984712 


1.00579 


18 


0.999997 


0.999816 


0.998332 


0.996139 


0.993357 


0.990678 


0.989654 


1.02348 


20 


0.999992 


0.999702 


0.997596 


0.994896 


0.991933 


0.989648 


0.997209 


1.04299 


22 


0.999983 


0.999551 


0.996756 


0.993664 


0.990808 


0.989309 


1.00689 


1.06337 


24 


0.99997 


0.999361 


0.995843 


0.992527 


0.990091 


0.989743 


1.01822 


1.08398 



TABLE III: Fitting (mn)" 1 \na {mn _ v)/2 to Eq. OJfor n = 11 
with 6 = 1 and cq fixed as Cg(ll) = 0.276016 given by Eq. 
^] Only data with m > mo = 100 are used in the fitting. 
Numbers in square brackets denote powers of 10. 
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-0.125966 


8.75546[-8] 


-8. 24421 [-6] 


0.000319302 


1.56222[-9] 


1 


-0.0561077 


-0.140372 


0.0319911 


-0.0345653 
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-0.0492583 


0.474695 


-4.7249 


23.4862 


0.999997 
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-0.0790665 
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-36.8528 


449.079 


0.999897 
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-0.133937 


3.79131 


-103.787 


1715.28 


0.999449 
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-0.207035 


6.24906 


-203.292 


3861.49 


0.998585 
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-0.293685 


9.05639 


-328.77 


6728.3 


0.997282 
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-0.391277 


12.162 


-480.3 


10421.5 


0.995726 
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14479.9 


0.993799 
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-0.609818 
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-833.317 
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0.991812 


10 


-0.727328 


22.305 


-1022.43 


24212.7 


0.989492 


11 


-0.850844 


25.9246 


-1235.21 


30014.7 


0.987267 


12 


-0.976389 


29.3893 


-1437.02 


35307.6 


0.984713 
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FIG. 1: (Color online) The original data of 
(mn)" 1 ln(a( mn _j,)/ 2 ) and the fitted curves for n = 11 
and v = 0, 1, . . . , 12. The curve for v — is on the bottom 
and that for v — 12 is on the top. The dashed horizontal line 
is Cq(11) = 0.276016 from exact expression Eq. 0] 
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FIG. 2: (Color online) The original data of 
(mn) -1 m(a( mn _„)/ 2 ) and the fitted curves for n = 12 
and v = 0, 2, ... , 24. The curve for v — is on the bottom 
and that for v — 24 is on the top. The dashed horizontal line 
is eg (12) = 0.279976 from exact expression Eq. 0] 



TABLE IV: Fitting (mn)" 1 In a (mn _ v)/2 to Eq. EJfor n = 12 
with b = 1/2 and c fixed as eg (12) = 0.279976 given by Eq. 
2] Only data with m > mo = 100 are used in the fitting. 
Numbers in square brackets denote powers of 10. 
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